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^ ' Abstract 

o ■ 

' We propose a new initial condition for the homogeneous and isotropic quantum cosmology, where 

, the source of the gravitational field is a conformally coupled scalar field, and the maximally symmetric 

fSJ ' hypersurfaces have positive curvature. After solving corresponding Wheeler-DeWitt equation, we 

' obtain exact solutions in both classical and quantum levels. We propose appropriate initial condition 

, for the wave packets which results in a complete classical and quantum correspondence. These 

^ ' wave packets closely follow the classical trajectories and peak on them. We also quantify this 
correspondence using de-Broglie Bohm interpretation of quantum mechanics. Using this proposal, 

^ . the quantum potential vanishes along the Bohmian paths and the classical and Bohmian trajectories 

^jq' coincide with each other. We show that the model contains singularities even at the quantum 

I— I, level. Therefore, the resulting wave packets closely follow the classical trajectories from big-bang to 
big-crunch. 

> : 

^ ■ 1 Introduction 
\0 '. 
<^ ■ 

• ' In recent years, the question of construction and interpretation of wave packets in quantum cosmology 

^ ' . . . . 

' and its connection with classical cosmology has attracted much attention. Moreover, some attempts 

o : 

^ ' have done to construct a theory of quantum gravity and find its relation with classical physics. Some 

' authors use the semiclassical approximations for the Wheeler-DeWitt (WDW) equation and consider 



the oscillatory or exponentially decaying solutions in the configuration space as classically allowed or 
forbidden regions, respectively. These regions are mainly determined by the initial conditions for the wave 
function of the universe. The no-boundary proposal of Hartle and Hawking [1] and the quantum creation 
of the universe from nothingness by Vilenkin are two popular proposals for the initial conditions. 

In quantum cosmology similar to the ordinary quantum mechanics, one is usually concerned with 
the construction of wave functions by the superposition of the energy eigenfunctions which would peak 
around the classical trajectories and follow them in configuration space, whenever such classical and 
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quantum correspondence is possible [31|31[S]. However, contrary to the ordinary quantum mechanics, an 
external time parameter is absent in the theory of quantum cosmology. Therefore, the initial conditions 
would have to be formulated with respect to intrinsic time parameter, which could be taken as the local 
scale factor for the three-geometry in the case of the hyperbolic WDW equation [6] . 

The construction of wave packets resulting from the solutions of the WDW equation is a common 
feature of some investigations in quantum cosmology [71 [HI [HI HH] ■ To be more specific, in Refs. [TUl [TTj 
the construction of wave packets in a Friedmann-Robertson- Walker (FRW) universe with a minimally 
coupled self interacting scalar field is presented and appropriate initial conditions are motivated. In 
Ref. [12j . the authors investigated (n + l)-dimensional quantum cosmology with varying speed of light 
and obtained exact solutions in both classical and quantum domains. The wave packets are constructed 
in such a way that they completely correspond to their unique classical counterparts. In Ref. [13] , a class 
of spherically symmetric Stephani cosmological models with a minimally coupled scalar field is studied in 
both classical and quantum levels [21 [T3] . The aim of these investigations has been to find wave packets 
whose probability distributions coincide with the classical paths obtained in classical cosmology. This 
means that the wave packet should be centered around the classical path and its crest should closely 
follow the classical trajectory. 

Here, we study a minisuperspace model describing a closed FRW universe with vanishing cosmological 
constant and containing a conformally coupled scalar field. The issue of the conformally coupled scalar 
field quantum cosmology has been in vogue for more than twenty years and basic results have been 
found by several authors [3l[4l[l6l[iail8l[l9l[2Ql|2ll[2ai23llM[25l|26l[^ In particular, Page [2^ has 
also studied this model and after solving the resulting WDW equation he obtained the solutions for the 
positive, negative and zero curvatures. For the case of positive curvature, the WDW equation casts into 
a oscillator-ghost-oscillator differential equation with well-known solutions. Moreover, Barbosa [26] has 
used one or two of these solutions to construct the wave packets and found the Bohmian trajectories via 
de-Broglie Bohm interpretation of quantum mechanics. Therefore, in contrast to the classical scenario, he 
found some non-singular solutions. In fact, the Bohmian trajectories highly depend on the wave function 
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of the system and various linear combinations of cigcnfunctions lead to different Bohmian trajectories. 
On the other hand, since the underlying WDW equation is second-order hyperbolic functional differential 
equation, we are free to choose the initial wave function and the initial slope of the wave function. But 
classically, we have a unique solution. Therefore, we encounter with a meaningful question: Is it possible 
to choose the initial condition in such a way that the resulting wave packet completely corresponds 
to its unique classical counterpart? In the previous investigations, using appropriate initial conditions 
(expansion coefficients), we could construct the wave packets which completely simulate their classical 
counterparts [TTl [121 US] . In particular, in Ref. [13], we obtained the wave packets so that the quantum 
potential vanishes along the Bohmian trajectories. In fact, the behavior of the quantum solution is strictly 
dependent to the initial conditions imposed on the wave function. Here, at the classical domain, the 
model is singular and we try to design wave packets in such a way that they follow as closely as possible 
the classical trajectories from big-bang to big-crunch. To achieve this purpose, we propose a specific 
relation between the expansion coefhcients or the initial wave function and its initial slope. Moreover, 
using WKB approximation we will show that the behavior of these wave packets is in agreement with 
the classical motion. 

The paper is organized as follows: in Sec. [2l we present the action of a homogeneous and isotropic 
quantum cosmology, where the source of the gravitational field is a conformally coupled scalar field, 
and the maximally symmetric hypersurfaces have positive curvature. In Sec. [31 we quantize the model 
and obtain the exact solutions of the corresponding WDW equation. We then construct the wave 
packets using the proposed initial condition. Moreover, using de-Broglie Bohm interpretation of quantum 
mechanics, we quantify the classical and quantum correspondence and justify the initial condition. In 
Sec. m we state our conclusions. 



Let us start from the Einstein-Hilbert action for the gravity plus a conformally coupled scalar field 



2 The model 




(1) 
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where is the four metric, g is its determinant, TZ is the scalar curvature and cj) is the scalar field. 
Units are chosen such that h = c — 1. Here, we consider a minisuperspace FRW model with the constant 
positive curvature and a homogeneous scalar field as 



ds2 ^ -.N^{t)dt^ +a^(t) 



l-r 



^ + {de'^ + sin^ 



(2) 



Now, we substitute Eq. in Eq. ([T]) and use the variable x = a^p0/%/2, where Ip is the Planck length 
(SttG = 3Zp). After discarding total time derivatives and integrating out the spatial degrees of freedom, 
we obtain the following action [201 ISS] 

(3) 



S = j dt 

This action results in the following Hamiltonian 

H 



aa^ ay , y 



a 



4a 4a 



^ y2 
a H 



(4) 



where Pa = — and P,(- = are the canonical momenta conjugate to a and x, respectively. These 

variables also satisfy the following Poisson brackets 

{a,x} = 0, {P„PJ = 0, 

(5) 

{a,Pj-l, {x,Px} = l. 
Using the above commutation relations and Hamiltonian equation (j4]), we find the equations of motion 

for the scale factor and the scalar field 

r d^{a,H} = -NPa/{2a), 



(6) 



Pa={Pa,H} =2N, 

X^{x.H} = NPj{2a), 

I Px^{Px.H} = ~2Nxla. 
By choosing the gauge N — a, we obtain the following parametric solutions for the system 

a(t) = Dsin(t), 

(7) 

x{t)^Dcos{t~eo), 

where D and 9q are constants. Also, we have used the zero energy condition for the super-Hamiltonian 
7Y = 0. It is obvious that above solutions represent Lissajous ellipsis which are singular in the present 
{t — 0) and in the future {t = tt). 



3 Quantum cosmology and wave packets 

Let us study the quantum cosmology aspects of the model presented above. The Hamiltonian can be 
obtained upon quantization procedure Pa — > ^-nd Px ~* Therefore, one arrives at the 

WDW equation describing the corresponding quantum cosmology [THl HOI US] 

n^ix,a) = l^~-^~X^ + aA^{x,a) = 0, (8) 



dx^ da?' 

where we have used a particular choice of factor ordering and for simplicity we have absorbed a factor 
of \/2 into variables. It is worth to mention that in the previous investigations [HI [El [E] , we have also 
encountered with the same form of WDW equation. There, we have constructed the wave packets in 
such a way that they corresponded to the classical cases where 0^ — 0. As we shall see, the usage of 
the real coefficients for the even terms and the imaginary coefficients for odd terms in the expansion, 
results in the symmetric solutions about a axis (^o — 0). Here, in order to study the non-symmetric 
cases (6'o ^ 0), we need to choose appropriate complex expansion coefficients for both even and odd 
terms. 

The WDW equation ([8|) is separable in the minisuperspace variables and solutions can be written as 

*„(x,a) = V„(x)V'«(a), (9) 



where 



e-^'/^ (10) 



In this expression, Hn{z) is the Hermitc polynomial and the orthonormality and completeness of the 
basis functions follow from those of the Hermite polynomials. The general wave packet which satisfies 
the WDW equation ([8]) can be written as 

*(x,a)= A^Mx)Ma)+i S„Vn(x)V'«(a). (11) 

n=even n=odd 

Since the potential in each direction is an even function, the eigenf unctions are separated in both even 
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and odd categories. Now, the initial wave function and its initial derivative take the form 

*(x,o)= AnMx)MQ), (12) 

even 

= * ^ B^MxWniO)- (13) 



da n 

n=odd 

Therefore, the coefficients A„ determine the initial wave function and the coefficients i?„ determine the 
initial derivative of the wave function. Since the underling WDW equation ([5]) is second-order hyperbolic 
functional differential equation, A^s and i?„'s are arbitrary and independent variables. On the other 
hand, if we are interested to construct the wave packets which simulate the classical behavior with known 
classical positions and momentums, all of these coefficients will not be independent. It is obvious that 
the presence of odd functions of a dose not have any effect on the form of the initial wave function but 
they are responsible for the slope of the wave function at a = 0, and vice versa for even functions. For 
studying the initial condition, let us consider the differential equation for the small values of the scale 
factor. Near a = 0, the WDW equation ([8]) takes the form 

This PDE is also separable in x and a variables, so we can write the solutions as 

ipn{x,a) =ipn{x)^n{a). (15) 
By substituting this expression in equation p4p . two ordinary differential equations can be obtained 



da? 



+ ^n^n(a) = 0, (16) 



+ x'Mx) = EnMx). (17) 



dx'^ 

where -En's are separation constants. These equations are Schrodinger-like equations with £^„s as their 
energy eigenvalues. The plane wave solutions are the exact solutions for Eq. (fTBl) 

£,n{a) = an cos (^y/E^ a^ + iPn sin (^^/E^ a^ , (18) 

where a„ and /3„ are arbitrary complex numbers. Equation (|17|) is the Schrodinger equation for the 
simple harmonic oscillator with the well known solutions (fTU)) . The general solution to equation can 



be written as 

V'(x,a) = J2 Kcos{y/Ka)tl;nix) + i K S'in{\/E^ia)Mx) ■ (19) 

n=even n=odd 

Note that, this sohition is vahd only for small a. Now, we can obtain the initial wave function and initial 
slope of the wave function 

V'(x,0) = ^nMx) (20) 

n— even 

^'(x,o) = i Y b:^,Mx), (21) 

n=odd 

where prime denotes the derivative with respect to the scale factor a. To have a complete description of 
the problem, we should specify both of these quantities. On the other hand, since we are interested to 
construct the wave packet with all classical properties, we need to assume a specific relationship between 

these coefficients. The prescription is that the coefficients have the same functional form [HI [121 [E] 

A* — C{n) for n even, 

(22) 

S* = C{n) for n odd, 
where C{n) is a function of n. In terms of AnS and i?„s we have 



An — I \^s C'{n) for n even. 



Bn = )v ^" C{n) for n odd. 



(23) 



Note that, we need to specify C(n) in such a way that the initial wave function has a desired classical 
description. We will see that this choice of coefficients results in a complete classical and quantum 
correspondence. Using equations pip and (|23p . we can explicitly write the form of the wave packet 



^" „iren^"(0) ^ ^ -^/^2"n! ^ 2nF„_,(0) ^(^^V4V2^ ^ 

(24) 

Although, in principle, we need to use infinite terms to construct the wave packet, for the studied cases 
a few terms are sufficient to get a reasonable accuracy |28| . Therefore, we use 50 terms in the above 
summation for all studied cases. 
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Figure 1: Left, the square of the wave packet a)P for C{n) = —^=e I""! 6*0 = tt/S and |C| = 4. 
Right, the contour plot of the same figure with the classical path superimposed as the thick solid line. 



m C'e 
Figure[T]shows the resulting wave packet for a particular choice of initial condition C{n) — — — j 



V2" n\ ' 

where C = ICIe"*^" ■ These coefficients are chosen such that the initial wave function consists of two 

well separated peaks. These two peaks correspond to classical initial {t = 0) and final (t = tt) values 

of X, respectively. As it can be seen from Fig. [1] the square of the wave packet is smooth and its crest 

follows the classical trajectory. In fact, we are free to choose any other appropriate initial condition. 

For example, figure [5] shows the resulting wave packet with different form of the expansion coefficients 

C{n) = — ^ ■ Figures [1] and [2] correspond to two different classical description with D = 4, 

V 2" ?i! 

6*0 = 7r/8 (Fig.[T]) and D = 5.7, = 7r/4 (Fig.©, respectively. 

Now, let us consider the situation using the ontological interpretation of quantum mechanics [29ll30j . 
This formalism allows us to compare quantum mechanical results with the classical ones. Moreover, 
since time is absent in quantum cosmology we can recover the notion of time using this approach. In 
de-Broglie Bohm interpretation of quantum mechanics, the wave function can be written as 

4' = i?e*^, (25) 

where R — R{x, a) and S — S{x, a) are real functions of the scale factor and the scalar field and satisfy 
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R dx^ R da? \dxj V . 

925 925. '2^dRdS_ _ ■2^dRdS_ 
g^2 Qg2 R dx dx R da da 



Figure 2: Left, the square of the wave packet \ip{x, a)P for C(n) = -^^^=e l^l^/"*, 6*0 = 7r/4 and \(\ — 5. 
Right, the contour plot of the same figure with the classical path superimposed as the thick solid line. 

the following equations 

Id^R Id'^R fdsV fdsV 2 2 

' - a' = 0, (26) 

0. (27) 

Because of the form of the wave packet (jlip . it is more appropriate to separate the real and imaginary 
parts of the wave packet 

'i^{x,a)=x{x,a)+iy{x,a), (28) 
where x, y are real functions of x and a 



-e 



2niJ„_i(0) ' ' " Tri/^VFd 



(29) 



^2ni7„_i(0) 7r^/*\/2"n! 



^=odd 

These variables are related to i? and S through equation (P5|) as 



(30) 



n— even 



i? = V^^ + y^, (31) 
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(I). ,32, 



arctan 

\x 

On the other hand, the Bohmian trajectories which determine the behavior of the scale factor and the 
scalar field, are governed by 

where the momenta correspond to the related Lagrangian ([3]) are ~ 2% and Pa = 2d, in the gauge 

N — a. Therefore, the equations of motion take the form 

dy dx 

2x = (35) 
dy dx 

2d = (36) 
x^ + 2/^ 

Using the explicit form of the wave packet (|29l30p , these differential equations can be solved numerically 
to find the time evolution of x ^nd a. In the left part of Fig. [3l we have shown the square of the wave 
packet |i/'(Xio)P for C{n) = — ^^=e~l''l^/^, 6q = 7r/8 and |C| = 7. In the right part of this figure, we 
have depicted x(0 ^"^^ o.{t) for classical (solid line) and Bohmian (dashed line) trajectories. In fact, the 
obtained Bohmian quantities versus time (x(0 ^^'^ '^i^)) coincide well with their classical counterparts. 
This shows the suppression of the quantum potential along the trajectory due to the coincidence between 
classical and Bohmian results [T^. 



For the case 0o 7^ 0, the momentum is not constant along the classical trajectory P = 2Dy sin^(t) + cos^(t — 6*0). 
In fact, we can also obtain some information about the classical momentum from the shape of the wave 
packet. For this purpose, we use WKB approximation. In semiclassical approximation, the square of 
the wave function up to the first order is related to the momentum as 

« ^- (37) 

This equation has a simple interpretation. The low momentum corresponds to the high probability 
density and vice-versa. Figure H] shows the square of the wave packet |i/'(X7^)P along the classical 
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Figure 3: Left, the square of the wave packet a)p for C{n) — ^ e 6*0 = 7r/8 and \(\ = 7. 

Right, plot of x{t) and a{t) for classical (solid line) and Bohmian (dashed line) trajectories. 
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t 

Figure 4: The inverse of classical momentum P^^ (dashed line) and the square of the wave packet 
|V'(Xi'j)P along the classical trajectory (solid line) for C{n) — — ^^=e~l^l = 7r/4 and |C| = 4 

versus t. 
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(Bohmian) trajectory which is parameterized with time and the inverse of the classical momentum P^^ 
versus time. As it can be seen from the figure, the height of the crest of the wave packet qualitatively 
shows the variation of the classical momentum along the trajectory. The low consistency between these 
two quantities is due to the approximate nature of equation (|37|) . 

4 Conclusions 

We have studied a closed homogeneous and isotropic quantum cosmology model in the presence of a 
conformally coupled scalar field. We have proposed a new initial condition which results in a complete 
classical and quantum correspondence. In fact, since the WDW equation is second-order differential 
equation, we are free to choose the initial wave function and the initial derivative of the wave function. 
This means that we are also free to choose the expansion coefficients. Since we are interested to have a 
consistency between classical and quantum solutions, we need to impose a particular relation between 
the coefficients. In this article, we proposed a particular relation between even and odd expansion 
coefficients. These coefficients determine the initial wave function and its initial derivative, respectively. 
In other words, this proposal defines a connection between position and momentum distributions with 
two properties at the same time. First, they correspond to their classical quantities and second, they 
respect the uncertainty relation. To quantify the classical and quantum correspondence, we have also 
studied the situation using de-Broglie Bohm interpretation of quantum mechanics. In fact, the Bohmian 
trajectories highly depend on the choice of the expansion coefficients and various linear combinations 
of eigenfunctions lead to different Bohmian trajectories. Therefore, although the inconsistency between 
classical and Bohmian solutions is natural in most cases, the quantum scenario is not always different 
from the classical scenario. In this paper, we have tried to construct the wave packets which peak around 
the classical trajectories and simulate their behavior. We showed that Bohmian positions and momenta 
coincide completely with their classical counterparts upon choosing arbitrary but appropriate initial 
conditions. Moreover, using WKB approximation, we qualitatively obtained the classical variation of 
momentum along the path without utilizing the classical equations of motion. 
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